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Abstract: Zarrineh River is located in the northwest of Iran, providing more than 40% of the total inflow 
into the Lake Urmia that is one of the largest saltwater lakes on the earth. Lake Urmia is a highly endangered 
ecosystem on the brink of desiccation. This paper studied the impacts of climate change on the streamflow 
of Zarrineh River. The streamflow was simulated and projected for the period 1992—2050 through seven 
CMIP5 (coupled model intercomparison project phase 5) data series (namely, BCC-CSM1-1, BNU-ESM, 
CSIRO-Mk3-6-0, GFDL-ESM2G, IPSL-CM5A-LR, MIROC-ESM and MIROC-ESM-CHEM) under 
RCP2.6 (RCP, representative concentration pathways) and RCP8.5. The model data series were statistically 
downscaled and bias corrected using an artificial neural network (ANN) technique and a Gamma based 
quantile mapping bias correction method. The best model (CSIRO-Mk3-6-0) was chosen by the TOPSIS 
(technique for order of preference by similarity to ideal solution) method from seven CMIP5 models based 
on statistical indices. For simulation of streamflow, a rainfall-runoff model, the hydrologiska byrans 
vattenavdelning (HBV-Light) model, was utilized. Results on hydro-climatological changes in Zarrineh River 
basin showed that the mean daily precipitation is expected to decrease from 0.94 and 0.96 mm in 2015 to 
0.65 and 0.68 mm in 2050 under RCP2.6 and RCP8.5, respectively. In the case of temperature, the numbers 
change from 12.33°C and 12.37°C in 2015 to 14.28°C and 14.32°C in 2050. Corresponding to these climate 
scenarios, this study projected a decrease of the annual streamflow of Zarrineh River by half from 2015 to 
2050 as the results of climatic changes will lead to a decrease in the annual streamflow of Zarrineh River 
from 59.49 m3/s in 2015 to 22.61 and 23.19 m3/s in 2050. The finding is of important meaning for water 
resources planning purposes, management programs and strategies of the Lake's endangered ecosystem. 
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1 Introduction 


Considerable changes in precipitation and temperature regimes in the past few decades, 
resulting in socio-economic and environmental issues, necessitates hydrological simulations for 
estimation of future changes in streamflow conditions. Global climate model (GCM), also 
known as general circulation model, using a mathematical model of the general circulation of a 
planetary atmosphere or ocean, can be taken into consideration to obtain climate change 
projections for climatic parameters such as temperature and precipitation (Liu et al., 2010, 2011; 
Yazdandoost et al., 2021). GCM simulations are based on the four global scenarios, RCP 
(representative concentration 
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pathways) 2.6, RCP4.5, RCP6.0 and RCP8.5, where RCP2.6 and RCP8.5 are the most optimistic 
and pessimistic scenarios, based on a possible range of radiative forcing values in 2100. Despite 
the increasing use of GCMs (global climate models), the main problem with them is their coarse 
resolution, which makes them unable to provide reliable information at the hydrological scales 
(Maraun et al., 2010). Therefore, using the downscaling approaches in relevant studies has become 
inevitable (Fowler et al., 2007; Teutschbein and Seibert, 2012; IPCC, 2013; McSweeney et al., 
2015; Bozkurt et al., 2018; Cai et al., 2018; Khan et al., 2018). 

Two major downscaling approaches, dynamical downscaling and statistical downscaling, are 
often used. Dynamical downscaling methods result in the higher-resolution regional climate models 
(RCMs), used to preserve the physical coherence between atmospheric and land surface variables 
(Anderson et al., 2011). However, the raw RCM (higher-resolution regional climate model) 
simulations may still be biased (Fowler et al., 2007). The statistical downscaling methods are based 
on statistical relationships between large-scale atmospheric variables (predictors) and local climate 
variables (predictands) (Do et al., 2011; Wilby et al., 2014). Compared to dynamic downscaling, 
statistical downscaling involves simple computational skills to downscale GCM simulations 
(Khalili et al., 2013). 

Artificial intelligence (AI), as a common method of statistical downscaling, is capable of 
analyzing long-series of hydrological data. In recent years, the use of AI methods such as artificial 
neural networks (ANNs) approaches is increasing. The success of this technique is primarily due 
to its ability to map highly nonlinear relations between inputs and outputs of the model (Coulibaly 
et al., 2005). Mendes and Marengo (2010) compared the performance of a temporal neural network 
model with an autocorrelation statistical downscaling model with emphasis on its ability to 
reproduce the observed climate variability and tendency for the period 1970-1999 in the Amazon 
region. The results showed that the neural network model remarkably outperforms the statistical 
models for the downscaling of daily precipitation variations. Skamarock et al. (2013) highlighted 
that ANN (artificial neural network) estimates of rainfall were more accurate than Weather 
Research and Forecasting model estimates, which applies dynamic downscaling approaches. 
Campozano et al. (2015) used the statistical downscaling model, the least squares support vector 
machines (LS-SVM) and ANNs to downscale precipitation data over the Paute River Basin in 
southern Ecuador. The results indicated that LS-SVM and ANNs performed better in downscaling 
of precipitation data. 

As development and implementation of adaptive management strategies in water systems 
requires the cognition of future changes in water resources, in this study, the aim was to simulate 
and project the effects of climate change on Zarrineh River's streamflow. In this regard, seven 
different climate model precipitation and temperature data series, were downscaled using an ANN 
technique and bias correction method. Then the bias corrected data were used in hydrological 
simulations of the streamflow. So, a toolbox was developed, which comprises the ANN toolbox of 
MATLAB (matrix laboratory), a program based on the bias correction method, a TOPSIS 
(technique for order of preference by similarity to ideal solution) program, and a rainfall-runoff 
model, namely the hydrologiska byrans vattenavdelning (HBV-light) model as the central engine. 


2 Materials and methods 


2.1 Study area 


Zarrineh River basin (35°40'-37°26'N, 45°46'—47°22’E), with an area of 1.18*10* km?, is located 
in the northwest part of Iran (Fig. 1). It provides more than 40% of the total inflow into Lake Urmia, 
which is on the brink of desiccation with a recorded annual water level reduction of about 40 cm 
over 1998-2018 (Ahmadaali et al., 2018). Lake Urmia, the largest lake in the country and the 
second largest saltwater lake on earth, has encountered a critical problem due to ineffective water 
resources management (comprising agricultural development and low irrigation water use 
efficiency, over allocation of the lake's inflows, irrigation projects and overuse of surface water and 
groundwater, etc.) and climate change (Fig. 2). 
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The mean annual precipitation and the annual mean temperature in the basin were reported to be 
about 292.00 mm and 13.26°C between 1992 and 2005, respectively. The potential daily 
evapotranspiration was about 4.17 mm. However, these reported numbers did not follow the 
specific trend during the period (Ahmadi et al., 2016; Ahmadaali et al., 2018; Alborzi et al., 2018; 
Yazdandoost et al., 2020). 

If a methodical comprehensive plan is built and implemented, reviving the shriveled lake is 
possible. As Zarrineh River feeds Urmia Lake, modelling the rainfall-runoff process in the river 
under the influence of future climate changes is crucial. The results are of important meaning for 
appropriate decision-making and water allocation in the basin. 
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Fig. 1 Location of the Zarrineh River basin in Iran (a) and the distributions of the meteorological and stream 
gauging stations in the study area (b) 


(a) 1993-1999 (b) 1999-2006 (c) 2006-2013 (d) 2013-2018 


Fig. 2 Urmia Lake shrinkage over 1993—2018 acquired by Land Remote-Sensing Satellite using Normalized 
Difference Water Index 


2.2 Data 


Zarrineh River basin, used as the case study of this work, was divided into twelve 0.5-degree 
cells. The cells were named from the top left. The study area lacks an adequate system to collect 
weather data. For this study, the observed daily mean temperature, provided by the Iranian 
Meteorological Organization (http://reports.irimo.ir/), were used as predictands at three observation 
sites, which are evenly distributed across the entire area of the basin (Table 1) and a bilinear 
interpolation technique was employed for developing gridded temperature data. Also, daily 
precipitation sums were provided from the global precipitation climatology center (GPCC), 
providing gridded gauge-analysis products, derived from quality-controlled station data. Here in 
the GPCC website (http://gpcc.dwd.de/), the gridded daily precipitation data is available in netCDF 
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(network common data form) format with a spatial resolution of 0.5°. Also, seven different CMIP5 
models for the most optimistic scenario, RCP2.6 and the most pessimistic scenario, RCP8.5, 
providing the predictor data, were gathered from the CERA (climate and environmental data 
retrieval and archive system) website (https://cera-www.dkrz.de/) (Table 2). For calibration of the 
HBV-Light (hydrologiska byrans vattenavdelning) model, measured streamflow values of Zarrineh 
River were used. More information of the four stream gauging stations has been provided in Table 
1. It should be noted that the calibration period of the study was 1992-2005, the evaluation was 
performed for the verification period 2006-2015, and the projection period was considered for 
2016-2050. 


Table 1 Location of the meteorological and stream gauging stations 


Station Location Altitude above sea level Period of 

Name Type Longitude Latitude (m) record 
Maragheh Meteorological 46°16'E 37°24'N 1477.7 1987-2019 
Saghez Meteorological 46°16'E 36°15'N 1522.8 1987-2019 
Takab Meteorological 47°06'E 36°24'N 1817.2 1992-2019 
Sarighmish Stream gauging 46°29'E 36°29'N 1372.0 1955-2019 
Gezkerpi Stream gauging 46°33'E 36°37'N 1350.0 1956-2019 
Miandoab Stream gauging 46°07'E 36°57'N 1290.0 1964-2019 
Nezamabad Stream gauging 46°03'E 37°03'N 1290.0 1993-2019 


Table 2 GCMs (global climate models) used in this study 


Resolution 

Model Complete name Atmospheric grid Ocean grid 
Longitude Latitude Longitude Latitude 
BCC-CSMI-1 ast ee 2.8125° 2.7906° noone SA 
BNU-ESM T a University-Earth 2.8125° 2.7906° 1.0000° ne 
cmomoso Comenta sented | grasa soe SET 
GFDL-ESM2G Pei Pa Dynamis 2.0000° 2.0225° Looe Pe 
IPSL-CMSA-LR ce e i theSth 3.7500° 1.8947° lon (ij) lat (ij) 
mocu Malice keh ags age aae OSE 
umoczsmcuem Anamo cheniay oip agas agoe aus OS 


Note: more information at https://portal.enes.org/data/enes-model-data/cmip5/resolution. lon (i, j) and lat (i, j) denote that the resolution 
cannot simply be read out. 


2.3 Artificial neural network 


ANNs are pattern recognition tools, which are able to find confirmable relationships between a set 
of inputs and the corresponding outputs (Chadwick et al., 2011; Patil et al., 2015). In 1957, 
Rosenblatt introduced the simplest form of feed-forward neural network, called perceptron. But as 
the inputs were fed directly to the output unit via the weighted connections in the method, it was 
proved that they could not be trained to distinguish different classes of patterns (Hung et al., 2009). 
A multilayer perceptron (MLP) is a feedforward ANN, comprising of at least three layers of nodes: 
an input layer, a hidden layer and an output layer. As the MLP finds non-linear relationships among 
the nodes, it overcomes many limitations of the single layer perceptron (Lippmann, 1987). Two 
important characteristics of the MLP are: its nonlinear processing elements which have a smooth 
nonlinearity; and their massive interconnectivity. The effectiveness of an ANN model is gauged by 
MSE (means quared error) and R?, which are derived from comparing the output data with the 
corresponding target output for the pattern. The best performance can be yielded by specifying the 
network architecture and by finding the optimal values of connection weights (Hung et al., 2009). 
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The network size depends on the number of nodes and layers. However, in the process of finding 
the optimal values for the connection weights, selection of a training algorithm is important. There 
are several training processes which differ in how the weights are obtained. In this study, the 
standard back propagation algorithm was used in ANN training based on its popularity. Detailed 
information about the proposed ANN approach is provide in Chadwick et al. (2011), TueVu et al. 
(2016), and Tran Anh and Taniguchi (2018). 

In the present work, to implement the ANN models, the ANN toolbox of MATLAB was used, 
and the transfer function of hidden layer was taken as tansigmoid while function of the output layer 
was a linear transfer function. The ANN models were trained on 70% dataset to predict precipitation 
and temperature for the twelve different 0.5-degree cells of the study area. The remaining 30% 
dataset was used for testing and validation of the models. In this study, it was desired to produce 
an output of temperature and precipitation and to establish the relationship between GCM 
simulations and observed temperature and rainfall data. ANNs were trained separately for 
temperature and precipitation, so each network produced only one output variable, corresponding 
to the target layer of the ANN. 


2.4 Bias correction 


Simulated precipitation and temperature by GCMs are often characterized by considerable biases 
due to their coarse spatial resolution (Gondim et al., 2018), limiting their direct application for 
basin-scale hydrological modeling. To address the systematic errors of GCMs, most of the climate 
studies use bias correction methods to bias correct the climate model outputs. This study describes 
how to use an efficient statistical bias correction method over the study area and improve the 
performance of the downscaled data, derived from the ANN models. The quantile mapping (QM) 
method (also known as distribution mapping method) was applied to correct precipitation and 
temperature data and is utilized to match the distribution function of the data, derived from ANN, 
with the observed distribution function. As it is used to adjust mean, standard deviation and 
quantiles, it preserves the extremes (ThemebI et al., 2012), but it also has its limitation due to the 
assumption that both meteorological variables follow the same distribution, which may cause new 
biases (Fang et al., 2015). The Gamma distribution (Thom, 1958) with shape parameter a and scale 
parameter J is often used for precipitation distribution (Piani et al., 2010): 


-x 


1 = 
=x"! ——— e’ (x20; a>0 0 1 
P ¥ e x20; œ >U, > 5 
Sap B°T(x) ( B ) ( ) 
where I(x) is the gamma function for the variable x. The QM method of precipitation is expressed 
mathematically in terms of the Gamma cumulative distribution function f, and its inverse fy | as: 


-1 
Porma = fy (i (Ponund | CANN, m 2 Pst) | A obs,m 2 Puc) 2 (2) 


where Peor,m,a 18 corrected precipitation; PANN,m,a is the precipitation derived from the ANN models 
on the day d of month m; aann,m and fann,m are the shape parameter and scale parameter from the 
ANN data on the month m; and Qobs, and fobs, are the shape parameter and scale parameter from 
the observational data on the month m. 

For temperature, the Gaussian (or normal) distribution (Cramer, 1999) with mean (u) and 
standard deviation (ø) is usually presumed to fit temperature best (Teutschbein and Seibert, 2012): 


1 =(x-ny 
= y% 20 u 
Ínu = * aie owe (3) 
where fy is the calculated Gaussian distribution and R is a real number. 
For temperature, the QM method is expressed mathematically in terms of the Gaussian 


r y É à à T = 
cumulative distribution function fy and its inverse fẹ as: 


-1 2 2 
Torma = Íy (fy tind | LLANN,m 2 O niin) | Hobs.m 2 Ocean) ? (4) 
where Teor,m,a is the corrected temperature; TANN,m,a is the temperature derived from ANN models on 
the day d of month m; HANN,m and 0” ANN.m are the mean and variance from the ANN data on the month 
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m; and [obs and 0 obs,m are the mean and variance from the observational data on the month m. 
2.5 TOPSIS method 


TOPSIS is a multi-criteria decision making method, which is based on the idea that the selected 
alternative should have the shortest geometric distance from the positive ideal solution and the 
longest distance from the negative ideal solution (Garcia-Cascales and Lamata, 2012). The method 
was first developed by Hwang and Yoon (1981). The procedure consists of the six following steps 
(Garcia-Cascales and Lamata, 2012; Jena et al., 2015): 

2.5.1 Establishing a performance matrix 

The decision-matrix (M) consists of m alternatives (A,,) and n criteria (C,), with the intersection of 
each alternative and criterion given as zj. 


G G a3, 

A Zi Zo a Diy 
M=A, | Zi Zg e Zm l (5) 

m Zm Zn 2 akti Z mn 


where i and j represent the number of row and column, respectively. 
2.5.2 Normalizing the decision-matrix 


In the classical TOPSIS approach, the normalized matrix can be obtained using the following 
transformation formula: 


N= (6) 


where Nj is a member of the normalized matrix. 
2.5.3 Calculating the weighted normalized decision matrix (V) as Equation 7, 
V= N; x W xn > (7) 


where these weights Wix, are calculated by entropy, AHP (Analytic Hierarchy Process) and direct 
assignation, etc. 


2.5.4 Determining the positive ideal solution (4*) and negative ideal solution (4°) 
A* = fvi, v} = {(maxv,, jE J\(minv,, je J’)}, i=1, 2, n m, (8) 
A = {v sV, y = {min v,, jE J)(maxv,, je J ^}, i=l, 2, ..., m, (9) 
where J is associated with benefit criteria; J’ is associated with cost criteria; and vy is the component 
of matrix V in the i line and j" column calculated from Equation 7. 


2.5.5 Calculating the separation measures 


The separation of each alternative from the A* (d;*) and the separation of each alternative from the 
A (di ) are given as follows: 


dj = lc DEPE? (10) 


Ba Oyo) 2 Pe eye, (11) 


2.5.6 Calculating the relative closeness (CL; ) to the ideal solution 
(so me (12) 
d; +d; 


Alternatives are ranked according to CL;. The higher CL; value, the higher the priority of the 
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i alternative (Garcia-Cascales and Lamata, 2012). 
2.6 Hydrological modeling 


HBV-Light, a conceptual, continuous and semi-distributed streamflow model developed by Sten 
Bergstrom in 1992 and 1995, was used to simulate daily streamflow values for the study area. HBV- 
Light includes some different routines and it can simulate snow, soil water, evaporation, 
groundwater and channel routing (Teutschbein and Seibert, 2012). HBV-Light requires daily 
temperature, precipitation and long-term average of potential evaporation values as driving 
variables. In this study, for the best GCM in the study area, chosen by TOPSIS method, 48 possible 
daily temperature and precipitation series in the twelve 0.5-degree cells, derived from the QM 
method based on the two scenarios of RCP2.6 and RCP8.5, were available to run the HBV-Light 
models. First, the HBV-Light models were calibrated and validated to measured streamflow values. 
For this purpose, the monthly averages of HBV-Light-simulated streamflow were analyzed and 
were compared to the observed values from 1992 to 2015. Then, these HBV-Light parameter sets 
were used to simulate streamflow with the corrected GCM simulations as driving conditions 
between 1992 and 2050. Figure 3 shows the flowchart of the study. 
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Fig. 3 Flowchart of the study. CMIPS5, climate model intercomparison project phase 5; GCMs, global climate 
models; RCPs, representative concentration pathways; TOPSIS, technique for order of preference by similarity to 
ideal solution. 


3 Results 


3.1 Statistical evaluation of the bias-corrected ANN simulations 


At this stage, we applied three statistical metrics to evaluate the performance of raw GCMs data 
over the study area and to assess the inter-model differences. The evaluation started by comparing 
the mean daily precipitation of GPCC data and individual model mean during the historical period 
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of 1992-2005. In addition, we used the mean absolute error (MAE) to compare the observed 
climate data with the corrected forecasted data and used the coefficient of determination (R°) to 
show the statistical relationship between the observation data set with the corrected data (Table 3). 
Results from Table 3 highlight the need for correcting GCMs outputs before any further analysis. 
We established the relationship between GCM simulations and observed temperature and rainfall 
data to implement the ANN models. In addition, we used the bias correction methods of 
precipitation and temperature to improve the performance of each ANN model. Precipitation and 
temperature correction methods were conducted on a daily basis from 1992 to 2050 (Table 4). 


Table 3 Statistical evaluation of the raw GCMs' simulations for daily temperature and precipitation data from 
1992 to 2005 in cell number 1 as a sample 


Temperature Precipitation 


Model i 
MAE (°C) R? Mean (°C) ar MAE (mm) R? Mean (mm) ae aa 
BCC-CSM1-1 6.95 0.68 15.12 1.21 0.02 0.60 
BNU-ESM 6.68 0.89 18.92 1.05 0.03 0.43 
CSIRO-Mk3-6-0 7.44 0.66 15.03 1.25 0.04 0.64 
GFDL-ESM2G 7.80 O71 18.53 13.26 1.36 0.01 0.74 0.80 
IPSL-CMSA-LR 6.35 0.77 922 1.30 0.02 0.73 
MIROC-ESM 6.57 0.90 19.38 1.25 0.01 0.64 
MIROC-ESM- 
CHIM 6.40 0.90 19.20 1.25 0.05 0.68 


Note: MAE, mean absolute error. 


Table 4 Statistical evaluation of the bias-corrected ANN (artificial neural network) simulations for daily 
temperature and precipitation data from 1992 to 2005 in cell number 1 as a sample 


Temperature Precipitation 
Model MAE (°C) R Mean (°C) nena is MAE (mm) R? Mean (mm) een 

BCC-CSM1-1 2.86 0.72 15,12 1.02 0.30 0.74 
BNU-ESM 2.77 0.93 18.92 0.87 0.30 0.74 
CSIRO-Mk3-6-0 2.36 0.69 15.03 1.09 0.32 0.77 

GFDL-ESM2G 2.85 0.76 18.53 13.26 1.21 0.30 0.77 0.80 
IPSL-CMS5A-LR 2.41 0.81 9.22 1.11 0.30 0.74 
MIROC-ESM 2.53 0.95 19.38 1.05 0.30 0.74 
MIROC-ESM-CHEM 2.54 0.94 19.20 1.06 0.30 0.74 


3.2 TOPSIS method 


The best CMIP5 model was chosen by TOPSIS method. Development of the TOPSIS algorithm was 
done based on the statistical evaluation of the bias-corrected ANN simulations from 1992 to 2005 
in the study area. As a non-limiting sample, Table 5 indicates that the highest priority of model 
selection is belonged to CSIRO-Mk3-6-0. Therefore, in the next step, this model was used to 
simulated Zarrineh River's streamflow under the RCP2.6 and RCP8.5 in the HBV-Light models. In 
this table, the bias-corrected ANN simulations for daily precipitation and temperature were 
evaluated statistically during 1992-2005. 

Table 6 shows the annual mean observed and corrected precipitation and temperature from 1992 
to 2005 by the bias-corrected ANN data, derived from CSIRO-Mk3-6-0 in Zarrineh River basin. 
Based on the table, a fairly good consistency between the observed and simulated data was revealed. 
Results indicate that the coefficient of determination between the annual mean observed and bias- 
corrected ANN data in the target period is 0.81 and 0.93 in the case of precipitation and temperature, 
respectively. 

Figure 4 indicates the observed and corrected mean precipitation and temperature from 1992 to 
2005 by bias-corrected ANN data for daily precipitation and temperature, derived from CSIRO- 
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Table 5 The result of TOPSIS (technique for order of preference by similarity to ideal solution) method for 
choosing the best CMIP5 (coupled model intercomparison project phase 5) model during 1992—2005 


Precipitation Temperature TOPSIS 

Mogel MAE (mm) R? Mean error (mm) MAE (°C) R? Mean error (°C) Cli* Rank 
BCC-CSM1-1 1.02 0.30 0.06 2.86 0.72 1.86 0.787 2 
BNU-ESM 0.87 0.30 0.06 2.87 0.93 5.66 0.106 5 
CSIRO-Mk3-6-0 1.09 0.32 0.03 2.36 0.69 1.77 0.971 1 
GFDL-ESM2G 1.21 0.30 0.03 2.85 0.76 5,27 0.289 4 
IPSL-CM5A-LR 1.11 0.30 0.06 2.41 0.81 4.04 0.450 3 
MIROC-ESM 1.05 0.30 0.06 2.53 0.95 6.12 0.028 7 
MIROC-ESM-CHEM 1.06 0.30 0.06 2.54 0.94 5.94 0.048 6 


Table 6 Observed and corrected mean annual precipitations and annual mean temperatures from 1992 to 2005 by 
bias-corrected ANN data, derived from CSIRO-Mk3-6-0 in Zarrineh River basin 


Time 
1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 
Pos (mm) 1.22 1.69 1.72 1.07 1.18 1.16 1.07 0.71 0.85 0.81 1.37 1.26 1.36 1.07 
Poor (mm) 1.30 1.34 1.29 1.19 1.13 1.08 1.01 0.96 0.96 1.01 1.09 1.13 1.12 1.00 
Tos (°C) 8.52 10.02 11.16 11.18 11.55 10.27 11.54 12.20 11.49 12.31 11.45 11.50 11.54 11.61 
Teor (°C) 9.76 11.10 12.12 12.56 12.42 12.20 12.46 12.86 12.95 12.96 12.87 12.74 12.58 12.45 


Note: Pops, observed precipitation; Peor, corrected precipitation; Tobs, observed temperature; Teor, corrected temperature. 
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Observed data —___— Corrected data 
84 (a) 


Precipitation (mm) 


Temperature (°C) 


Number of cells 


Fig. 4 The observed (points) and corrected (box plots) mean precipitation (a) and mean temperature (b) from 
1992 to 2005 by bias-corrected ANN data for daily data, derived from CSIRO-Mk3-6-0 


Mk3-6-0. In this figure, the cells were named from the top left. 

In the validation period of 2006-2015, bias-corrected ANN data for precipitation and 
temperature data of CSIRO-Mk3-6-0 projects, based on RCP2.6 and RCP8.5, were assessed based 
on the observed data. The coefficient of determination in the validation period for bias-corrected 
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ANN data for precipitation was 0.50 and 0.51 in the case of RCP2.6 and RCP8.5, respectively. 
These numbers change to 0.70 and 0.74 when it comes to temperature. As Figures 5 indicate, results 
achieved an acceptable agreement between the annual mean observed and bias-corrected ANN data. 

Bias-corrected ANN data for precipitation data of CSIRO-Mk3-6-0 reveals that the daily mean 
precipitation is expected to decrease from 0.94 and 0.96 mm in 2015 to 0.65 and 0.68 mm in 2050 
based on RCP2.6 and RCP8.5, respectively. In the case of temperature, the numbers change from 
12.33°C and 12.37°C in 2015 to 14.28°C and 14.32°C in 2050, respectively. This data was next 
used to simulate Zarrineh River's streamflow in the HBV-Light models. 
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Fig. 5 Assessment of bias corrected ANN simulations of monthly mean precipitation (a) and temperature (b) 
derived from CSIRO-Mk3-6-0 in the validation period (2006-2015) 


3.3 Evaluation of the hydrological models 


The chosen CMIP5 model (CSIRO-Mk3-6-0), operated under RCP2.6 and RCP8.5, was used to 
simulate Zarrineh River's streamflow in the HBV-Light models from 1992 to 2050. For the 
evaluation of the models, statistical indexes such as the coefficient of determination were used. 
Also, monthly mean streamflow in hydrological simulations for the catchment were compared with 
the observed data between 1992 and 2015. In general, it can be said that the agreement between the 
observed and simulated values both in the calibration period of 1992-2005 and in the validation 
period of 2006-2015 is acceptable. The observed monthly streamflow in 2015 and corresponding 
simulated streamflow in HBV-Light models by the corrected data, derived from CSIRO-Mk3-6-0, 
operated under RCP2.6 and RCP8.5, were summarized in Tables 7 and 8. Generally, all simulations 
were in good agreement with the observations and there was no significant difference between these 
results and the observed values. 

Upon assuring that the models simulate the stream flow of Zarrineh River in the plausible limits, 
future projections of the river flow were achieved. Future projections for the period of 2016-2050 
of streamflow for the study site, simulated by bias-corrected ANN data for daily precipitation and 
temperature, derived from CSIRO-Mk3-6-0, operated under RCP2.6 and RCP8.5, were generally 
in agreement in terms of the general trend (Fig. 6). A more detailed look at the figure reveals that 
the annual streamflow of Zarrineh River is expected to decrease from 59.49 m?/s in 2015 to 22.61 
and 23.19 m?/s in 2050 based on RCP2.6 and RCP8.5, respectively. 
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Table 7 Performance of the HBV-Light (hydrologiska byrans vattenavdelning) models simulated by bias- 
corrected ANN data for daily bias-corrected precipitation and temperature, derived from CSIRO-Mk3-6-0, operated 
under RCP2.6 and RCP8.5 in 2015 


R? Sum of monthly mean streamflow (m°/s) Monthly mean streamflow (m°/s) 
RCP2.6 RCP8.5 Observed RCP2.6 RCP8.5 Observed RCP2.6 RCP8.5 
0.63 0.63 713.84 703.58 714.34 59.49 58.63 59.53 


Table 8 Monthly mean streamflow of the HBV-Light models simulated by bias-corrected ANN data for daily 
bias-corrected precipitation and temperature, derived from CSIRO-Mk3-6-0, operated under RCP2.6 and RCP8.5 
in 2015 


Monthly mean streamflow (m*/s) 


Month Simulation 
Observation 
RCP2.6 RCP8.5 
Jan 42.74 32.40 34.70 
Feb 211.78 38.10 37.75 
Mar 202.61 219.21 219.81 
Apr 89.31 75.21 75.01 
May 45.76 54.20 53.72 
Jun 40.69 45.50 45.42 
Jul 27.30 42.89 42.89 
Aug 16.68 40.82 40.83 
Sep 9.89 39.18 39.36 
Oct 7.83 37.22 38.74 
Nov 8.19 33.77 36.23 
Dec 11.06 45.08 49.88 
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Fig. 6 Annual mean streamflow of the Zarrineh River in the HBV-Light (hydrologiska byrans vattenavdelning) 
models simulated by bias-corrected ANN data, derived from CSIRO-Mk3-6-0, operated under RCP2.6 and RCP8.5 
from 2015 to 2050 


4 Discussion 


Climate change is expected to alter the magnitude and timing of runoff and, consequently, have 
significant implications for existing water resources systems as well as for planning and 
management of future water resources (Teutschbein and Seibert, 2012). Zarrineh River basin faces 
high pressure on water due to ineffective water resources management, as well as climate change 
(Abbaspour and Nazaridoust, 2007; Ahmadaali et al., 2018; Alborzi et al., 2018). So, a better 
understanding of the impact of climate change on water resources in the basin is necessary 
(Ahmadzadeh et al., 2011; Delju et al., 2013; Nikbakht et al., 2013; Gholampour et al., 2015; Torabi 
Haghighi and Klove, 2017; Ghale et al., 2018). The present study aimed to simulate and project the 
streamflow of Zarrineh River from 1992 to 2050. For this purpose, a toolbox was developed, 
comprising the ANN method for downscaling the large-scale data, a bias correction method for 
generating unbiased climate datasets, a TOPSIS program for choosing the best CMIP5 model and 
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a rainfall-runoff model, namely HBV-Light for investigating the effects of climate change on water 
resources. The results indicated that the daily mean precipitation is expected to decrease from 0.94 
and 0.96 mm in 2015 to 0.65 and 0.68 mm in 2050 based on RCP2.6 and RCP8.5, respectively. In 
the case of temperature, the numbers change from 12.33°C and 12.37°C in 2015 to 14.28°C and 
14.32°C in 2050, respectively. In addition, the annual streamflow of Zarrineh River is expected to 
decrease by half from 2015 to 2050. 

Results from this study, along with previous studies, showed that the combination of the ANN 
method and bias correction method can improve the performance of the forecasted data. As the 
results showed that the performances of bias correction methods are not robust in all of the months. 
To evaluate the robustness of the performances of the methods in different hydrological seasons, 
the streamflow can be divided into different periods according to the hydrographs. As statistical 
downscaling methods rely heavily on the assumption that currently observed relationships will 
carry into the future, using two statistical downscaling methods simultaneously may result in 
similar consequences for both RCP2.6 and RCP8.5. 

The proposed approach can be further strengthened and extended from several aspects. Future 
research will, hopefully, combine and apply the ANN method and different bias correction methods 
to other study areas as the accuracy of the method depends on several spatial and temporal 
conditions. A major extension may be considering land cover-use change, where any change in land 
cover may also result in a decrease in streamflow while sensitivity analysis is an approach to handle 
this type of challenge. It is undeniable that other potential topics can be covered in future studies 
in Zarrineh River basin to quantify the impacts of different future scenarios and to reduce the 
negative effects of climate change and human activities on local water resources. 


5 Conclusions 


Considerable changes in precipitation and temperature regimes in the past few decades have 
resulted in socio-economic and environmental problems. Accordingly, it is required to estimate 
future changes, being used in hydrological simulations of streamflow. In spite of the increasing use 
of GCM simulations in hydro-climatic studies, their direct application is contestable due to the 
biases. This paper introduces a new approach for forecasting rainfall and temperature, using an 
ANN technique. To improve its performance, a bias correction method may be used. So, a toolbox 
was developed, comprising the ANN method, a bias correction and a TOPSIS program, and a 
rainfall-runoff model, namely the HBV-light model, where the effects of climate change on water 
resources were investigated. The toolbox was utilized for the case of Zarrineh River basin. Future 
projections of streamflow revealed that the annual streamflow is expected to decrease by half from 
2015 to 2050. 
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